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Abstract 

Conformation of single polyelectrolytes in tetravalent salt solutions is investigated under the 
framework of a coarse-grained model, using Langevin dynamics simulations. The chain size, studied 
by the radius of gyration, shows three different variational behaviors with salt concentration, 
depending on the chain stiffness. According to the size variations, polyelectrolytes of fixed chain 
length are classified into three categories: (1) flexible chain, for which the variation shows a curve 
similar to a tilted letter L; (2) semiflexible chain, whose curve looks resemble of the letter U; (3) 
rigid chain, for which the curve is a straight line. The worm-like chain model with persistence 
length predicted by the Odijk-Skolnick-Fixman theory is found to be able to qualitatively describe 
the end-to-end distance at low salt concentration, not only for semiflexible and rigid chains but also 
for flexible chain. In a low-salt region, a flexible polyelectrolyte extends more significantly than 
a semiflexible chain, in reference of the size of their uncharged counterparts, and in a high-salt 
region, regardless of chain stiffness, a chain attains a dimension comparable to that of its neutral 
polymer. The chain stiffness influences both the local and the global chain structures. A flexible 
chain exhibits a zigzagged local structure in the presence of salt ions and the condensed structure is 
a disordered, random globule. A semiflexible chain is locally smooth, and the condensed structure 
is orderly packed, taking a form such as hairpin or toroid. Moreover, the chain stiffness can also 
affect the nature of the coil-globule transition. The transition is occurred in a discrete manner 
for semiflexible chain, whereas in a continuous way for flexible chain. This discrete feature is 
happened not only at low salt concentration when a semiflexible chain is collapsed, but also at 
high salt concentration when the collapsed chain is reexpanded. At the end, the effects of chain 
stiffness and salt concentration on the conformation of single polyelectrolytes are summarized in a 
schematic state diagram. 

PACS numbers: 
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I. INTRODUCTION 



There is a resurgent interest in studying the properties of polyelectrolytes in multiva- 
lent salt solutions because such system reveals many fascinating phenomena. One vital 
example is the DNA condensation, in which DNA, a negatively-charged polyelectrolyte, un- 
dergoes a dramatic condensation from an extended structure to a compact, highly-ordered 
structure while multivalent salt is added to the solution [H [2l |3l H]. The condensation in- 
duced by multivalent salt is not a privilege reserved only for DNA but a common feature 
of polyelectrolytes (HJ El E] ■ Usually trivalent salts or charged molecules of higher valence 
are demanded to induce the condensation [21 E]- Experiments have shown that the mor- 
phology of the condensate depends strongly on the chain bending rigidity which balances 
the attractive collapsing force at a solvent quality. Flexible polymers generally collapse to 
disordered globules, and stiff ones collapse to ordered structures such as toroid or folded 
chain [21 [71 E]. The influence of chain stiffness on chain morphology has been investigated 
by simulations, including the condensation of polyelectrolyte induced by multivalent salt [TU] 
or by increasing Coulomb strength parameter [TT] and the collapse of neutral polymer driven 
by temperature decrease [121 [13]. The results confirmed that chain rigidity plays a decisive 
role in determination of the condensed structure. 

A phenomenon which is less well-known appears when an excess of multivalent salt is 
presented in the solution: The salt-induced condensation redissolves into the solution and 
the system returns to a homogeneous phase [51 151 [E] . This phenomenon is called "reentrant 
condensation" [15] and has been discovered since seventy-five years ago in the study of 
the coacervation of colloids by oppositely charged ions [16] . It occurs for flexible chains, 
for example, polystyrene sulphonate 0, and also for stiff chains, for example, DNA [6]. 
The persistence length in the latter example is two order of magnitude larger than in the 
previous. As the condensation takes place, the bare chain charge is almost neutralized by the 
condensed counterions dissociated from the multivalent salt [171 [IB] . These counterions not 
only screens out the long-range Coulomb repulsion but also creates a short-range attraction 
due to correlated fluctuation, resulting in the condensation [191 120] . On the other hand, for 
the redissolution of polyelectrolytes at high salt concentration, no comprehensive explanation 
is well established. Recently, two theories were proposed to explain the phenomenon but 
verification is still ongoing. The first theory suggested that the redissolution is a consequence 
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of charge inversion occurring when polyelectrolyte binds so many counterions that its net 
charge alters sign [ini EI], whereas the second one argued that it has an origin of the 
screening of the short-range attraction The second theory, furthermore, predicted 

that Bjerrum association plays an important role on the organization of polyelectrolyte-ion 
complexes [21], and therefore, charge reversal is not always happened. 

The folding transition of a polymer between an elongated state and a compact state is 
called "coil-globule transition" [25]. Theorists predicted that for semifiexible chains, such 
transition is first order, using mean-field theory |2S1 1251 12Z]- The prediction has been con- 
firmed experimentally [2Sl [2H]- Despite continuous for the whole ensemble, the coil-globule 
transition is discrete at the level of single chain for polymers with sufficient stiffness |29] . In 
this case, the transition is similar to the disordered-ordered transition between a gas-like and 
a crystal-like state. On the other hand, for single polymers with small stiffness, the transi- 
tion is continuous, resembling to the disordered-disordered transition between a gas-like and 
a liquid-like state. This topic has been recently studied in simulations [IHl [HI [121 [13] ■ The 
results showed that the discrete transition can be induced by adding condensing agent [10] 
or by decreasing temperature [121 [H], while chains are stiff. However, the discussions were 
mainly focused on the chain condensation occurring in low temperature region or in low 
salt regime. Since an excess of multivalent salt can lead condensed chains reentering into 
the solution, it becomes very relevant to know whether the reentering transition behaves 
in a discrete manner or not. To our knowledge, there is no simulation nor experiment, 
hitherto, discussing the conformational transition on the chain reentrance occurring at high 
salt concentrations. To get insight of it, we investigate in this paper the behavior of single 
polyelectrolytes with addition of multivalent salt by means of computer simulation. The salt 
covers a broad range of concentration so that single chains show, in turn, the condensation 
and the reentrance transitions in a microscopic way. Our study, therefore, offers a good 
opportunity to investigate thoroughly the coil-globule transition happened at low and at 
high salt concentrations. 

There have been many simulations devoting to the study of the behavior of polyelec- 
trolytes in salt-free solutions [IIl[30l[ai[321[33l[32[35l[36l[37l|38l[Mlll^ Only recently, the 
study with addition of salt becomes numerically feasible because of the progress of comput- 
ing power. However, the salt valence studied in the literatures was usually monovalent or 
divalent [ITl WI\ HH [S]. For the study with high salt valence, the focus was mainly on the 
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condensation of polyelectrolyte and the salt concentration was not high enough to investi- 
gate the chain reentrance into a solution [TUl HHl SSI HZ] • One indirect way to discuss DNA 
condensation and redissolution is to measure the effective interaction between two immobile 
rigid polyelectrolytes in the environment of multivalent salt [391 H8] . Recently, progress has 
been made in a series of study [IHl EDI El] : By carefully choosing the simulation parameters, 
we were able to investigate salt-induced condensation and redissolution thoroughly in a di- 
rect way. We found that the dimension of flexible polyelectrolytes undergoes subsequently 
two continuous transitions upon addition of multivalent salt, the collapsing transition and 
the swelling transition, which are respectively the microscopic representation of the conden- 
sation and the redissolution of polyelectrolyte. We demonstrated that the excluded volume 
of ion plays an imperative role on the chain properties and should be inevitably incorpo- 
rated in a theoretical analysis |l9l [50]. The study of the potential of mean force showed 
that like-charge attraction between chains happens only when salt concentration is inter- 



mediate around the equivalence point (see Sec. Ill A for definition) and the size ratio of ion 
to monomer lies within a window around unity ^U\. Moreover, it was shown that only salt 
valence greater than two can induce apparently the two structural transitions [51j, which 
is consistent with experimental observations [H |6l [171 E2]- Nonetheless, the effect of chain 
stiffness was not discussed in these studies. Since semiflexible polyelectrolytes exhibit several 
intriguing phenomena, not being seen in flexible polymers, we focus in this work how chain 
stiffness influences chain morphology. 

Toroid structure is presumed to be the ground state of a compacted DNA [21 [S3]. How- 
ever, folded-rod structure and a number of other morphology, such as racquet, are also often 
observed in experiments [5H [551 [56] . The long lifetime of some of these structures indicates 
that DNA has complicated energy landscape. Numerical studies showed that the energy 
levels of toroid and rodlike structure are comparable [121 [25] . These two structures are sur- 
rounded by many metastable states and hindered by high energy barriers. Debate continues 
until now in deciding which structure is the real ground state of a condensed DNA and under 
what condition it applies. One of the important factors which affects the chain morphology 
is the chain length. It was observed that longer chains display higher appearance percent- 
age for toroids than for rods [21 [571 [SH]. Also, the degree of ion condensation on a chain 
was found to increase with chain length [SHI [SD] • Recently, chain length effect on the state 
diagram of a single chain has been discussed [SI]. Limited to today's computational power. 
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time period able to be studied by molecular simulations is very short. It is, thus, difficult to 
investigate properly the evolution of chain morphology, because of the long lifetime of some 
of the chain structures. Therefore, many authors found that the formation of a toroidal or a 
rod-like structure pertains to the initial chain configuration used in simulations pTl [T2| l35] . 
In this article, we study the morphology of a single polyelectrolyte with fixed chain length. 
For the chain stiffness in which the morphological evolution is slow, several simulations are 
performed, starting with different initial configurations, to make sampling more complete. 
By this way, we calculate static properties. The rest of the article is organized as follows. 
The simulation model and the method used in this study are described in Sec. [TTj The results 



and discussions are given in Sec. Ill The topics include the radius of gyration (Sec. Ill A) 



the degree of chain swelling (Sec. IIIB), the end-to-end distance (Sec. Ill C[ ), the probability 



density distribution (Sec. HIE), and the dynamics of chain conformation (Sec. IIIF). Sim 



ulation snapshots are presented in Sec. [Ill D We summarize the conformation of a single 



polyelectrolyte in a state diagram in Sec. Ill G and give our conclusions in Sec. IV 



II. MODEL AND SIMULATION METHOD 

Our system consists of a single polyelectrolyte modeled by a negatively charged bead- 
spring chain and spherical counterions dissociated from the chain. There are 48 beads 
(monomers) on the chain. Each of the beads carries a — e charge and each of the coun- 
terions carries a +e charge. Tetravalent salt is added into the system. It is dissociated 
into tetravalent cations (counterions) and monovalent anions (colons), modeled by charged 
spheres. Solvent molecules are not incorporated exphcitly in the simulations. We suppose 
them forming a medium of constant dielectric constant e,.. The collision between particles 
and solvent molecules is considered, and the effect is described by Langevin equation. The 
system is placed in a cubic box and periodic boundary condition is applied. 

Four types of interaction are involved in the simulations. The first one is the excluded- 
volume interaction, which is applied for all the particles, including monomers, counterions 
and colons. It is modeled by a purely repulsive Lennard- Jones (LJ) potential 

, , / [{a/rr - [cr/rf] + e,. for r < 2^6^ 
Uex{r) = \ (1) 
I for r > 2^l^a 

where Elj is the interaction strength and a is the diameter of a particle. The second 



interaction is the Coulomb interaction Ucoui{f) — ZiZjc^ / {Aireoerr) where Zi and Zj are 
respectively the valences of particle i and j, eo is the permittivity in vacuum, is the 
dielectric constant of the medium. It can be re-expressed as 

Ucoui[r) = 2 

r 

where ks is the Boltzmann constant, T is the temperature, and = / {AireoerkBT) is 
the Bjerrum length, which denotes the distance between two unit charges at which the 
Coulomb energy is equal to the thermal energy ksT. The third interaction is the bond 
connectivity interaction. Two adjacent monomers are connected by a spring, modeled by a 
finitely extensible nonlinear elastic potential 

Ubondib) = - 2^6&max In f 1 - — J (3) 

where h is the length of a bond, h^ax is the maximum bond extention, and k\) is the spring 
constant. The forth interaction is a harmonic angle potential, which gives chain stiffness, 

Uan,le{Q)^k,{e-Q^f (4) 

where Q is the angle between two consecutive bonds, is the equilibrium angle, and ka is 
the force constant. 

We perform Langevin dynamics simulations in this study. The equation of motion of a 
particle % is described by the stochastic differential equation 

mjfj = - mi-iiVi + ^i{t) (5) 

where is the mass of particle i, mi7i denotes the friction coefficient and -f^i is a random 
force modeling the collisions by solvent molecules. 77*1 (i) has zero mean over time and satisfies 
the fiuctuation-dissipation theorem: 

(^i(^) • = 6kBTmai5ij5{t - t') (6) 

where 6ij and 6{t — t') are the Kronecker delta and the Dirac delta function, respectively. 
An advantage of using Langevin dynamics is that the temperature control is naturally in- 
corporated. 

In this work, we assume that all the particles have identical mass m, diameter a, and 
LJ interaction strength Elj- We set = 3a and control temperature at /c^T = 1.2eLj- 
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Coulomb interaction is computed using the technique of Ewald sum. kb is chosen to be 
5.8333kBT/(T'^ and bmax to be 2a. Under this setup, the average bond length is l.la [311 
HH]. We study the effect of chain stiffness by altering the spring constant ka from to 
100kBT/ra(P . The equilibrium angle 6q is set to be vr. The size of the simulation box 
is 53.133cr and the concentration of tetravalent salt ranges from 0.0 to 0.01024cr~^. All 
the particles are subjected to a friction force proportional to the particle velocity. We 
set the damping constant 7 to 15. Or^^ to mimic roughly an aqueous environment where 
T = a^Jm/ {kBT) is the time unit. Simulation time step 5t is chosen equal to O.OOSr. The 
initial chain configuration begins with an extended structure. An equilibrium phase takes 
3 X 10^ to 5 X 10^ time steps and a production run 10^ to 3 x 10*^ time steps, depending 
on chain stiffness and salt concentration. We found that in the region of intermediate salt 
concentration, the transition between different chain structures is slow. In order to sample 
more completely different structures, five to ten independent runs are performed starting 
with different initial configurations. The simulations were run using LAMMPS package [7B] . 
For the case of the highest salt concentration, the system contains 7776 charged particles. 



III. RESULTS AND DISCUSSIONS 



A. Radius of gyration 

The size of a polyelectrolyte can be characterized by the radius of gyration, Rg. It is 
calculated by taking the square root of the formula, 

1 ^ 
1=1 

where is the number of monomers on the chain, is the position vector of monomer i, 
and f*cm is the center of mass of the chain. We studied the effects of salt concentration and 

1 /2 

chain stiffness on the root-mean-square radius of gyration, (-Rg) • The results are shown 
in Figjl] 

We observed that the variation of (-Rg) behaves in three distinct manners. The first 

1 /2 

behavior applies to the chain with small fca, such as ka = O.OkBT/racf. (-R^) firstly 
decreases with increasing Cs, the concentration of the tetravalent salt, until the total charge 
of the tetravalent cations is equal to the negative charge carried on the chain backbone. We 
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call this particular salt concentration the equivalence point. In this work, the equivalence 
point occurs at = C* = 8 x lO^^a^"^ because the monomer concentration is fixed at 
Cm = 3.2 X lO^^cr^^. (-Rg) then increases slowly while Cs surpasses the equivalence point. 
Therefore, the chain undergoes two structural transitions upon addition of tetravalent salt: 
a shrinking transition and a followed swelling transition. The increasing rate of (i?^)^''^ in 
the chain swelling region is apparently lower than the decreasing rate in the chain shrinking 

1 /2 

region. As a result, the curve of {Rg) looks similar to a tilted 'L' in the semilog plot. 
Notice that at high salt concentration, the chain is not swollen back to its original size in 
the absence of salt. These phenomena have been recently studied in detail by molecular 
dynamics simulations |l9l [50] . 

The second behavior applies to the chain with intermediate ka, for instance, ka = 
8.0kBT/ra(P . At the beginning of addition of salt, the chain size decreases gently. The 
gentle decrease turns to become a sharp decrease once Cg is increased beyond some critical 

1 /2 

value Cc, and quickly, {R^) reaches a small constant value. If Cs is further increased 

1 /2 

over a second critical value Cd, an abrupt increase in {Rg) takes place. The equivalence 
point C* is situated in the region between Cc and Cd- In the high-salt region Cs > Cd, 
the chain reattains roughly its original size before the condensation occurred. These two 
sudden changes suggest that the chain morphology is altered, from an extended structure to 
a condensed structure, and vice versa, in a discontinuous way. The chain possibly commits 

1 /2 

a first-order conformational transition [121 [22] • In the semilog plot, the shape of the (R^) 
curve looks like the capital letter 'U'. A relevant phenomenon has been observed in the 
experiments of DNA precipitation by polyamines where DNA percentage presented in the 
supernatant shows a similar U-shaped variation as a function of spermidine concentration 
[6l [62]. Suppose that the occurrence of the chain precipitation is directly related to the 
chain size in the solution. This variation will then refiect how the dimension of a semifiexi- 
ble DNA varies with salt concentration. Our results are in agreement with the experimental 
observation. 

The third behavior is observed for the chain with large ka, for example, ka = 

1 /2 

100kBT/ra(P. The infiuence of the adding salts upon the chain size is not significant. (-Rg) 
stays essentially at a constant value, and in the figure, the curve is roughly a straight line. 

Based upon the variation of chain size against Cg, a polyelectrolyte of fixed chain length 
can be classified into three categories according to the chain stiffness: 
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1 /2 

(1) flexible chain, characterized by a tilted-L-shaped (-Rg) curve; 

(2) semiflexible chain, characterized by a U-shaped (i?^)^^^ curve; 

1/2 

(3) rigid chain, characterized by an almost straight (-Rg) curve. 

One way to justify if a chain belongs to the category rigid chain or not is to compare the 
bare persistence length with the chain length. If the bare persistence length is larger than 
the chain length, the chain is considered to be rigid. In our study, the bare persistence length 
ip^o is related to the force constant ka of the harmonic angle potential by ip^ = 2ka{b) / {kBT) 
where (6) is the average bond length. Therefore, ^pfl becomes larger than the chain length 
while ka > 24kBT/rad'^. This simple criterion works reasonably well because the chain 
shows a rigid rod structure at all Cg for ka = SOksT/rad'^ and lOOfc^T/rad^ but not for the 
cases with smaller ka, as having been shown in Figjl] 

In the simulations, we observed that, no matter how stiff the chain is, the polyelectrolyte 
is in its most stretched state while no salt is added into the solution. We found that at 

1 /2 

a given salt concentration, a stiffer chain displays a larger (-R^) , and hence, occupies a 
larger space. The results obtained here also demonstrate that the chain stiffness plays an 
imperative role in the determination of the type of the coil-globule transition to be continuous 
or discrete, both in the condensation and in the reexpansion of a polyelectrolyte. 

B. Degree of swelling 

In order to understand the degree of swelling of a polyelectrolyte in salt solutions, we 
calculate the size of the uncharged counterpart of the polyelectrolyte in a salt-free solution 
and use it as the reference of comparison. The root-mean-square radius of gyration of the 
uncharged chain, (i?^ is equal to 5.1(1), 12.1(2), and 14.97(6) for the cases ka = 
0.0, 8.0 and lOO.OkBT/rad"^, respectively. We define the degree of chain swelling by a = 
{{RlY^^ I (^g,„y^^) - 1 and present the resuhs in Fig. |2| 

We observed that in the low-salt region, a is positive. It means that the charged chain 
occupies a larger space than the uncharged chain, apparently resulting from the electrostatic 
repulsion between monomers on the chain. The value of a for the flexible chain (k,a = 
O.OkBT/rad"^) is much larger than the ones for the semiflexible chain {ka = S.OkBT/rad'^) 
and for the rigid chain [ka = lOO.OfcsT/raci^). Therefore, the major contribution to the 
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dimension of a semiflexible or a rigid polyelectrolyte comes from the intrinsic chain stiffness, 
whereas the size of a flexible chain is mainly determined by the Coulomb repulsion on the 
chain backbone. In the mid-salt region around the equivalence point, the degree of swelling 
is negative for the flexible chain and for the semiflexible chain. The chains are, hence, in 
collapsed states. In this region, electrostatics dominates the chain stiffness. The maximum 
shrinkage, in reference of the size of the neutral chain, can be as large as 50%. For the 
rigid chain, the collapsing transition does not happen, and thus, a is close to zero. In the 
high-salt region, we saw that the value of a approaches to zero for the semiflexible chain. 
The polyelectrolyte hence occupies a space volume similar to its uncharged counterpart. 
For the flexible chain, a is increased with Cg but still takes a finite negative value. In our 
previous study [50j, flexible polyelectrolytes in the region of much higher salt concentration 
have been investigated, and the results showed that the chain attains a size similar to the 
neutral polymers. Therefore, a is expected to tend toward zero if was increased to a 
much higher value in the current study. 

In Fig. |3| we present the variation of a as a function of at different Cg- At small 
Cs, a is positive and monotonically decreases to zero with increasing ka- Opposite trend of 
behavior is observed at large C^, for instance, Cg = 0.01024o"~'^, where a is negative and 
monotonically increases to zero. In the mid-salt region, a displays a V-shaped curve and the 
maximum shrinkage is happened around ka = bksT /ra(P . Thus, by tuning ka and varying 
Csi one is able to control the degree of chain swelling to design desired functions in real 
applications. 

The results obtained here tell us that in a salt-free solution, the larger the value of fca, the 
less significant the size difference between a charged chain and its uncharged counterpart 
will be. Since at large Cg a chain occupies a space similar to the one of a neutral polymer, a 
semiflexible chain will reattain roughly its original size at zero Cg by addition of many salts 
into the solution, but a flexible chain will not reattain it. 

C. End-to-end distance 

Conventionally, the end-to-end distance Re is used to characterize the size of a polymer. In 
Fig. [4] we present how the root-mean-square end-to-end distance, (Rl)^^'^ = {{vn — n)^)^''^, 
varies with salt concentration for different chain stiffness, where fi and vn are the position 
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vectors of the two ends of a chain. The behavior of Re resembles the behavior oi Rg-. a 
flexible chain displays a tilted L-shaped Re curve, a semiflexible one displays a U-shaped Re 
curve, and a rigid one displays a straight Re curve. Nonetheless, the value of Re does not 
precisely reflect the real size of a chain. For example. Re for a semiflexible chain becomes 
smaller than that for a flexible chain in the mid-salt region, although a semiflexible chain has 
a larger size (cf. Fig. [T|. This inconsistency results from the formation of ordered structures 
which gives a small Re-, such as hairpin or ring-like structure. Therefore, the quantity Re is 
not always suitable to characterize the dimension of a chain in a salt solution, particularly 
when the chain stiffness is intermediate. 

The mean-square end-to-end distance of a chain has been derived theoretically under the 
framework of the worm-like chain (WLC) model and reads as 

{Rl) = 2l,L - 2il (1 - expi-L/Q) (8) 

where L is the contour length and ip is the persistence length [63]. For a polyelectrolyte 
chain, ip is the sum of two contributions: the first one is the bare persistence length ip^, 
coming from the intrinsic chain stiffness, and the second one is the electrostatic persistence 
length ie, due from the electrostatic interaction. Odijk [64J and, independently, Skolnick 
and Fixman (OSF) have deduced a formula for the electrostatic persistence length in 
the limit of stiff chain and low ionic strength. This formula reads as 



ffOSF 



(9) 

where ^ is the Manning parameter, defined as the ratio of over the effective charge 
distance on the chain, and k is the inverse Debye length. In order to make a comparison of 
the simulation results with the theory, we adopted the Manning-Oosawa (MO) condensation 
theory [191 ES] to estimate ^. In our case, ^ is renormalized to a value 1/4, due to the 
condensation of tetravalent counterions, while Cs > C* = Cm/ 4:. For small salt concentration 
Cs < C*, all the tetravalent counterions are expected to condense on the chain, but the 
number of these ions is not sufficient to renormalize ^ to 1/4. At this moment, ^ decreases 
linearly from 1 to 1/4 with Cs, or equivalently, ^ = 1 — 3Cs/Cm- Plugging ^ and = 
AnXBiCni + 4C, + A^Cs) into Eq. g, we obtained £f^^, and then, added it to ip^ = 
2ka{b) /{kBT) to get the persistence length ip. Finally, using Eq. (8), we calculated (Rl)^^"^ 
and the results for different chain stiffness are plotted in dashed curves in Fig. |4j We observed 
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that the WLC model, combining with the OSF and the MO theories, gives a quahtative 
prediction for {R'^)^^'^ at low salt concentration for semiflexible and stiff chains. Surprisingly, 
it also gives a not-bad description of (-R^)^^^ for the flexible chain with zero ka and even 
captures the decreasing behavior up to Cg = C*, although the OSF and the MO theories are 
not expected to work well neither for flexible chains nor for a solution with multivalent salt. 
The theoretical value is found to be smaller than the value obtained from the simulations. 
This deviation can be attributed to the effect of the excluded volume interaction, which 
is neglected in the derivation of the WLC model and should give a positive contribution 
to the chain size if it is taken into account. In the mid-salt regime, the theory fails to 
predict (Rl)^^'^ for semiflexible chains. In this region, the chain is condensed by tetravalent 
counterions. The small chain size suggests that the persistence length is significantly reduced, 
and as a consequence, the electrostatic persistence length takes a negative value. The 
negative electrostatic persistence length has been experimentally observed [67]. At high 
salt concentration, the WLC model again describes well the end-to-end distance except for 
the case with zero ka, at this moment, i^^^ is almost zero and the persistence length is 
dominated by the bare persistence length. 

D. Snapshots of simulation 

Before processing a more quantitative analysis, we present in this section snapshots of our 
simulations to give readers a vital illustration how the chain looks like in different regions 
of salt concentration and chain stiffness. 

Fig. [5] shows snapshots of a flexible chain at low, middle, and high salt concentrations. 
The morphology of the chain is elongated at low salt concentration, and displays large-angle 
turns near the places where the counterions condense, as shown in Fig. [sj^a). In the mid-salt 
region around the equivalence point, the chain collapses and exhibits a randomly-arranged 
compact structure (see Fig.|5](b)). At high salt concentration (Fig. [5|^c)), the chain swells and 
acquires a zigzagged, extended structure, but not as elongated as at low salt concentration. 

Fig. [6] shows snapshots of a semiflexible chain at three salt concentrations. The shape of 
the chain is smoother than that of the flexible chain. At low salt concentration (Fig. [6]^a)), 
the chain is elongated and not zigzagged. At middle salt concentration, the chain is con- 
densed, and noticeably, forms ordered structure, in contrast to the disordered globule formed 
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by a flexible chain. Typical ordered structures are shown in Fig. [6[b). They are, from left 
to right, hairpin, racquet, and toroid. These typical structures have been observed in the 
experiments of DNA condensation [71 \5E\. Unlike the case for flexible chains, a 

condensed counterion can hardly induce a large-angle turn on the chain backbone. If it is 
the case, the turn happens to be the folded end of a hairpin structure and counterions are 
condensed between the two branches of the hairpin. The whole structure looks similar to 
peas in a pod. At high salt concentration, the chain displays again an elongated structure 
(see Fig. gc)). 

Fig. [7] shows the snapshot of a rigid chain. The chain exhibits a stretched structure no 
matter how many salts are added into the solution. The chain is stiff enough to resist from 
bending or twisting. 

From these snapshots, we found that a compact chain structure is formed when condensed 
tetravalent counterions can overcome the chain stiffness to bridge between monomers and 
have enough number to tightly bind the whole chain structure. We observed that, regardless 
of the chain stiffness, the number of tetravalent counterions condensed on a chain increases 
with Cs- When Cg is smaller than the equivalence point C*, almost all the tetravalent 
counterions are condensed on the chain. Only when Cg > C*, non-condensed tetravalent 
counterions can appear in the bulk solution. Moreover, the number of condensed tetravalent 
counterions surpasses the number needed to neutralize the bare chain charge when Cg > C*, 
and as a consequence, the chain is locally overcharged [77]. For flexible chains, there are 
mainly two effects to drive chain reentrance at high salt concentration. The flrst one is the 
Coulomb repulsion between condensed tetravalent counterions. Since the chain is surcharged 
by the condensed tetravalent counterions, the net repulsion between these ions manifests 
its effect in the swelling of the chain size. The second one is the presence of the non- 
condensed tetravalent counterions in the bulk solution. These non-condensed counterions 
interact with the chain, and can dynamically replace previously condensed counterions and 
become newly condensed counterions. During this process, the chain morphology oscillates 
from a loose globule to an expanded structure and vice versa. Increasing salt concentration 
increases the concentration of non-condensed tetravalent counterions in the bulk, which 
makes easier the replacement of tetravalent counterions to take place. The chain hence 
more frequently stays in an expanded state and the averaged chain size increases. For 
semiflexible chains, there is, in addition, a third effect taking part in the reentrant transition 
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— the chain stiffness. The chain stiffness always gives an effect against the formation of a 
compact structure. Together with the instabihty induced by the presence of non-condensed 
tetravalent counterions, a chain can dramatically leave a condensed state to an elongated 
structure at high salt concentration. 



E. Probability density distribution 

In the previous sections, the conformation of a polyelectrolyte was studied by averaged 
quantities, such as the root-mean-square Rg and the root-mean-square Re. In fact, averaged 
quantities cannot reflect precisely the chain shape since chain morphology varies with time 
during simulations. In order to understand the details of chain morphology, we investigate 
in this section the probability density distribution (PDD) of conformational quantities. 

Two quantities are studied: The first one is the asphericity A and the second is the radius 
of gyration Rg. The asphericity A measures the deformation of chain morphology from a 
spherical geometry and is defined by [68] 

Ai, A2, and A3 are the three eigenvalues of the gyration tensor of the chain calculated by 

1 ^ 

Taf3 = ^(^i - fcm)airi - fcm)f3 (H) 
1=1 

where the subfixes a and (3 denote the three Cartesian components x, y and z. The value 
of A ranges between and 1. It is equal to for a perfect sphere, 0.25 for a perfect ring, 
and 1 for a straight line. For a coil chain, (A) is 0.431 obtained from simulation |69j . 

In Fig. [sj parts (a) and (b), we present PDD of asphericity, p{A), and PDD of radius of 
gyration, p{Rg), for the flexible chain with null ka. We observed that p{A) shows a peak 
near A = 0.8 in the salt-free solution, indicating that the chain favors an extended structure. 
The peak gradually moves toward a small A as Cg is increased, and the width of the peak 
becomes broader. The width then decreases while Cg approaches to the equivalence point 
C* (Case VI in the figure). At C*, p{A) shows the most pronounced peak at A ~ 0.1, 
and therefore, the favored morphology is a compact sphere-like structure. If Cs is further 
increased, the peak becomes broaden and the position of the peak shifts not much toward a 
larger value. In this salt region, the chain behaves like a random coil and alternates between 
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a globule and an extended structure. For p{Rg), similar behavior was observed. The peak 
of p{Rg) moves left to a small Rg and then moves slightly right as Cs is increased. Also, the 
most pronounced peak appears at the turning point Cg = C*. 

Figs.[9]^a) and (b) shows, respectively, p{A) and p{Rg) for the semiflexible chain with ka = 
SkBT/rad"^. In the low-salt region, p{A) displays a single peak at A ~ 0.9, and therefore, 
the chain favors to exhibit a rod-like structure. If the salt concentration is increased to 
some intermediate value, a second peak appears at A ~ 0.25. It indicates the formation of a 
ring-like (toroid) structure. The second peaks first increases in its height with Cs and then 
decreases. It disappears while Cg surpasses some critical value, and the chain reattains a 
rod-like structure. More information can be obtained if we study p{Rg) at the same moment. 
For example, at Cs = 2.67 x lO^^a^^ (Case III in the figure), p{Rg) shows a doubly peaked 
distribution and the second peak is roughly situated at the middle position between the 
first peak and Rg = 0. We know from p{A) that the chain exhibits a rod-like structure at 
this salt concentration. The second peak in p{Rg), hence, corresponds to hairpin structure 
because hairpin is a two-folded rod and its radius of gyration is equal to half of that of 
the unfolded rod. Consequently, the chain morphology alternates between unfolded rod and 
hairpin structure. At higher salt concentration, the first peak in p{Rg) disappears and a 
new peak appears at Rg ^ 5cr. This new peak is attributed to toroidal structure because it 
corresponds to the peak of p{A) at A ^ 0.25, which is the typical value for a ring structure. 
Therefore, in this salt region coexist two condensed structures, hairpin and toroid. We noted 
that probability density for A or Rg is nearly zero between the two peaks when Cs is not large, 
which indirectly indicates that the structural transition between hairpin and toroid rarely 
happens during a run. If Cs is large, for instance, Cs = 6.4 x 10~'^a~^ , a finite value appears 
between the two peaks, indicating that the two structures can frequently transit between 
each other. We mention that these results were obtained by doing statistics combined 
from several independent runs starting with different initial configurations to circumvent 
the problem of slow transition. At Cs = 0.00128cr^'^ (Case X in the figure), p{Rg) shows 
three peaks, indicating a coexistence between the two condensed and the extended-chain 
structures. One thing is worth to be noticed: In this middle salt region, the peak in p{Rg), 
which corresponds to a toroidal structure, shifts toward a large value as Cs is increased. 
It suggests that the toroidal size grows up with salt concentration. This phenomenon is 
in agreement with the observation of Conwell et al. [70j. In the experiment, they found 
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that monovalent or divalent salt causes an increase in diameter of a condensed DNA toroid. 
Finally, if Cg goes even higher, p{Rg) displays only one peak at Rg ~ 13a, and hence, the 
favored shape is an unfolded chain. 

We emphasize that the PDDs for the semiflexible chain shows different behavior, com- 
pared with the ones for the flexible chain. In the middle salt region, they are multiply peaked 
and the peaks appear or disappear in a sudden way. Moreover, each peak corresponds to 
an ordered structure such as toroid or hairpin. On the other hand, the PDDs for the flex- 
ible chain display only single peaks and the peak position moves continuously in the space 
with salt concentration. As a consequence, the coil-globule transition for the semiflexible 
chain occurs in a discrete manner, whereas in a continuous way for the flexible chain. Both 
the discrete and the continuous transition have been demonstrated in simulations by other 
authors [10], [121 [121 US]. However, the subjects were focused on the chain collapsing from a 
coil state to a globule state, happened at low salt concentration or at low temperature. Our 
study presented here goes further and gives a thorough description of how the PDDs varies 
with salt concentration, not only in the low salt region but also in the high salt region. The 
sudden appearance and disappearance of peaks show that the discrete nature of transition 
takes place also at high salt concentration when a semiflexible chain redissolves from an or- 
dered globule state to a coil state. To our knowledge, this is the first time such observation 
at high Cg is reported. 

We present in Fig. [Io|(a) and Fig. [Io|(b) p{A) and p{Rg), respectively, for a rigid chain 
with ka = lOOksT/racP. p{A) shows a sharp peak near A = 1 and p{Rg) a sharp peak at Rg 
around 15a for all the salt concentration investigated. The chain exhibits a fully extended 
structure. The results can be verified by estimating the chain length from the radius of 
gyration. It is known that the length £ of a rod is related to the radius of gyration by 
R^ = i'^/12. Consequently, Rg = 15a yields £ equal to 51.96(t, which is in good agreement 
with the length of a fully extended chain since our chain has 47 bonds and the average bond 
length is l.lcr. 

We remark that the results presented here were obtained from long simulations. We have 
calculated the correlation time for the quantities A and Rg at each salt concentration. The 
results shows that our simulations contains hundreds to thousands independent data for the 



cases of flexible chain and rigid chain. Therefore, we are confident that Fig. M and Fig. 10 



report correct equilibrium populations. For the case of semiflexible chain, particularly in 
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the mid-salt region, the system shows rich morphologies and the transition between different 
morphologies is slow. We have preformed several independent runs to sample different 
trajectories in order to improve the statistics. The analysis of the correlation time shows 
that at least, thirty independent samples were cumulated in the simulations. The results 
reported in Fig. |9j hence, may not be very accurate. However, the PDDs show consistent 
trend of variation against C^, and also in Fig. [T| no remarkably large fluctuation is appeared 

1 /2 

in the (-R^) curve. Therefore, these results still give an appreciable description of the 
equilibrium populations. 



F. Dynamics of conformation 

We have shown that a polyelectrolyte can be condensed by multivalent salt. Particularly, 
when the chain is semiflexible, the condensed chain displays an ordered structure, such as 
hairpin or toroid. In order to get insight of the phenomena, we take the semiflexible chain of 
ka = SkBT/rad"^ as an example, and study in this section how the conformational quantities, 
A, Rg and Re, evolve with time. 

We focus at first the conformational dynamics at = 8 x lO^^a^^. In Fig. [9| we 
have demonstrated that the PDDs of the semiflexible chain are doubly peaked at this salt 
concentration. The chain stays in a hairpin state or a toroid state and the transition between 
the two states is not frequent. For this case, several independent runs were performed 
starting with different initial configurations. In order to get insight of the dynamics, we 



show firstly in Fig. 11 the time evolution of A, Rg and Re for one of the independent runs, 



in which the chain exhibits hairpin structure. 



The asphericity A fluctuates steadily around 0.9 as shown in Fig. 11 a). At the same 
moment, the radius of gyration fluctuates around Rg = 6.5a and Re fluctuates around a 
small value (see Fig. [lT|^b)). Therefore, the chain has its two ends fluctuating near the open 
end of the hairpin structure. Readers can see from the figures that this structure is lasted 
for a duration of at least 2 x 10^ simulation steps in this run. Please be aware that 2 x 10^ 
simulation steps is a long time period in a typical molecular simulation. In our case, it is 
approximately equal to a duration of 2/xs, provided that a is 2.4A, T is 300K and m IS 
200Da. 



We then show in Fig. 12 the dynamics of these quantities when the chain exhibits toroidal 
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structure. The asphericity A fluctuates around 0.25 except at four places, fluctuating around 



0.7 (cf. Fig. 12 ^a)). At these four places, the chain displays a twisted ring structure, similar 
to the shape of the digit '8'. Since the two ends of a toroidal chain is generally not very 
close to each other. Re fluctuates around a large value (see Fig. [T2|b)), compared to the 
value for hairpin structure (Fig. [lT|b)). Resembling the previous run, the toroid structure is 
continued for nearly 2 x 10^ simulation steps. The morphological transition between toroid 
and hairpin is, hence, not frequent at this salt concentration. 

The transition between different ordered structures becomes observable in simulations if 
the salt concentration is high. For example, at Cg = 6.4 x 10~'^a~^ (shown in Fig. [l3|a)), 
the asphericity fluctuates around two values, aX A = 0.25 and at A = 0.85, which indicates 
the transition between ring-like and rod-like structures. Moreover, Fig. [I3|(b) shows that 
the radius of gyration associated to the rod-like structure is larger than that associated 
to the ring-like structure, but the end-to-end distance is smaller. It enables us to identify 
this rod-like structure to be a two-folded chain [i.e. hairpin), which was also confirmed 
from the simulation snapshots. The second example is chosen at Cs = 1.28 x lO^^a^^ (see 



Fig. 14). The asphericity shows similar fluctuation around the two values. However, the 
value, A ~ 0.85, which corresponds to rod-like structure, is now associated to two behaviors 
of Rg, fluctuating at Rg ~ 6.5a and Rg ~ 13cr. Combining the information obtained from 
Re, we assert that the first fluctuation of Rg corresponds to a two-folded chain structure and 
the second one corresponds to an unfolded- chain structure. Therefore, the chain transits 
between three states, which are unfolded chain state, hairpin state, and toroid state. 

Recently, dynamic exchange between toroid and folded-chain structures in a salt-free 
solution of semiflexible polyelectrolyte has been numerically demonstrated by tuning the 
Coulomb strength parameter |TT]. The work showed a back-and-forth flipping between the 
two structures, similar to our observations. Our study, moreover, showed that transitions 
between different chain states can be enhanced by addition of multivalent salt beyond the 
equivalence point C*. The mechanism is understood as follows. Below C*, most of the 
multivalent counterions are condensed on the chain. Roughly while Cg is larger than C*, 
multivalent counterions start to exist in the bulk solution. Therefore, increasing salt concen- 
tration increases the number of the multivalent counterions surrounding but not condensed 
on the chain. Due to the presence of these surrounding counterions, the chance for a chain to 
escape from its current condensed state is increased. The surrounding multivalent counteri- 
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ons strongly interact with the chain, and can induce local change of the chain morphology by 
the replacing process: some of the condensed ions on the chain are released into the solution 
and replaced with the surrounding ions. Occasionally, the chain forms a newly condensed 
structure. 



G. Schematic State diagram 

The conformation of our single polyelectrolytes, classified by the three categories: flexible 
chain, semiflexible chain, and rigid chain, is summarized in a schematic state diagram shown 



in Fig. 15 



Flexible chain exhibits an extended structure in a salt-free solution, and continuously 
collapse to a compact, spherical structure as Cg is increased to the equivalence point. Further 
increasing Cs gradually leads chain discondensation. At very high salt concentration, the 
chain behaves similar to a random coil. 

Semiflexible chain shows different conformations as multivalent salt is added into the 
solution. At very low salt concentration, the chain exhibits an elongated structure. If Cs 
is increased beyond the first critical value Cc, the chain bends over upon itself and forms 
two- folded rod structure (hairpin). This structure is not so stable that the chain morphology 
alters between a hairpin and an extended structure. Further increasing Cg firstly stabilizes 
the hairpin structure, and eventually, leads to the appearance of toroidal structure. At this 
moment, hairpin and toroid structures coexist in the solution. We found that the rate of 
transition between hairpin and toroid increases with the salt concentration. Therefore, the 
transition becomes observable during limited time period of a molecular simulation while 
Cg is large enough. If Cg is not large enough, the two structures can be observed only in 
independent runs, since our model system contains solely one chain. The formation of these 
ordered structures has been proposed to be a consequence of a nucleation-growth process 
derived from the thermal fluctuations in chain morphology [7T1 EH]. At higher salt 
concentration, the chain can even transit between three morphological states: the unfolded- 
rod state, the hairpin state, and the toroid state. Finally, if Cg surpasses the second critical 
value Cd, the chain can no longer be collapsed, and consequently, reattains a unfolded, 
rod-like structure. 

Rigid chain, just as its name denotes, is rigid, no matter how many salts are added in 
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the solution, and always shows a rodlike, fully stretched structure in simulations. 

The state diagram presented here is in accordance with the results obtained by other 
authors, although the employed models are different. For example, Ou and Muthukumar 
investigated the conformation of single polyelectrolytes in a salt-free solution and the chains 
were collapsed by increasing Coulomb strength parameter F [TT]. They observed that a 
semiflexible chain stays, in turn, in a coil state, a folded-rod state, and a toroid state, as F 
is increased (cf. Fig. 6 of Ref. [H]). The result is consistent with what we observed in the 
region where chains are collapsed. For the second example, we mention the works done by 
Ivanov and his coworkers [121 EI]- In their works, the effects of chain stiffness and chain 
length on the conformation of neutral polymers were studied. They found that finite chain 
length can make unsharpened the coil-to-toroid transition driven by temperature reduction, 
and several intermediate states, such as folded-rod and disk-like structures, were observed in 
the transition region. The proposed state diagram is in analogy of our results in the low-salt 
region. We remark that chain length is an important factor in determination of the chain 
morphology. Experiments have shown that short DNA molecules are more easily condensed 
to folded rods than toroids [21 [571 EH]- This phenomenon gives an explanation of why the 
hairpin structure appears frequently in our state diagram. It is because the modeled chain 
is not long in our study. Nonetheless, the schematic state diagram obtained here has gone 
beyond the works done by other authors. We investigated the salt-induced conformation 
of polyelectrolytes thoroughly in a direct way since this topic has attracted much attention 
in scientific community [7^ [75] , and explored the situations from a salt-free condition to a 
very high salt concentration so that chain condensation and redissolution are both occurred. 
The findings give deep insight of the phenomena of DNA condensation at molecular level 
and are relevant for a broad range of salt-induced complexation phenomena. 



IV. CONCLUSION 



We have performed Langevin dynamics simulations to study single polyelectrolytes in 
tetravalent salt solutions under the framework of a coarse-grained model. The effect of salt 
concentration and the role of chain stiffness on the properties of chain conformation have 
been discussed. In order to understand thoroughly the reentrant condensation, a broad 
range of salt concentration has been investigated so that the single polyelectrolytes display. 
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in turn, collapsing transition and swelling transition upon addition of salt. 

We have observed that chain morphology crucially depends on chain stiffness. At a fixed 
chain length, the radius of gyration Rg of a polyelectrolyte varies in three different ways 
with the salt concentration Cg- Based upon these variations, we can classify polyelectrolytes 
into three categories. The first category is flexible chain for which (i?^)^^^ shows a tilted 
L-shaped curve as a function of Cs in the semilog plot. The second one is semiflexible chain 

1 /2 

which is characterized by an U-shaped (-Rg) curve. And the third one is rigid chain for 
which (R^) is roughly a constant, no matter how many salts are added in the solution. 
The end-to-end distance shows the similar variation as the radius of gyration. However, its 
value does not always reflect the real chain size, particularly when the chain is condensed to 
an orderly packed structure such as hairpin. Moreover, we have found that the WLC model 
with persistence length obtained by the OSF theory qualitatively predicts the end-to-end 
distance at low and at high salt concentrations for semiflexible and rigid chains. Surprisingly, 
these theories also give a not-bad description of the decreasing behavior of (Rl)^^^ against 
Cg for flexible chain up to Cg — C*, although they are not expected to work well neither for 
flexible chains nor for a solution with multivalent salt. 

We have found that in a salt-free or vary low-salt solution, the degree of swelling of a 
polyelectrolyte decreases with increasing the chain stiffness, in reference of the size of the 
uncharged counterpart of the chain. On the other hand, in a high-salt region, it tends 
toward zero, regardless of the chain stiffness. As a consequence, at high Cg, a semiflexible 
chain reattains a dimension close to the chain size in the absence of salt, but not does 
a flexible one. We have discussed the reasons for the occurrence of chain reentrance or 
swelling at high salt concentration. For flexible chains, two effects can possibly drive this 
phenomenon: (1) the net Coulomb repulsion between condensed multivalent counterions 
which locally overcharges the chain, (2) the instability derived from the presence of non- 
condensed multivalent counterions in a bulk solution. For semiflexible chains, the chain 
stiffness plays, in addition, the imperative role on reexpanding the chain. 

The chain stiffness has been witnessed to influence both the local and the global structures 
of a polyelectrolyte. A flexible chain shows a zigzagged local structure in the presence of 
multivalent salt and the condensed structure is disordered. In contrast, a semiflexible chain 
locally shows smooth morphology. Large-angle turn on the chain backbone can only be 
observed when the chain is condensed to a hairpin structure. Moreover, the global structure 
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of a condensed chain is orderly packed. Hairpin and toroid are the two favored structures. 

The chain stiffness affects also the nature of the structural transition. By studying the 
probability density distributions of asphericity A and of Rg, we have demonstrated that 
a semiflexible chain undergoes a coil-globule transition in a discrete manner between an 
extended structure and a condensed, ordered structure, whereas a flexible chain undergoes 
a continuous transition. This discrete nature is happened not only in the low-salt region 
when chain condensation occurs, but also in the high-salt region when chain reenters into 
the solution. 

We have studied the dynamics of the conformational quantities of a semiflexible chain. 
Upon addition of salt beyond the first critical concentration Cc, the chain is folded to a hair- 
pin structure. However, it is a dynamic equilibrium process and the chain morphology alters 
between a hairpin and a unfolded-chain structure. Further increasing Cg firstly stabilizes 
the hairpin structure and then triggers the birth of toroid structure. Transition between 
hairpin and toroid becomes more frequent as getting higher. When Cg approaches the 
second critical concentration C^, these two structures become unstable and the chain can 
stay again in a unfolded state, coexisting with the two condensed structures. At the end, we 
summarized the conformation of individual polyelectrolytes in a state diagram. Although 
the modeled chain is not very long, our results agree with experiments and other numerical 
works in many aspects in the chain condensation region. The work presented here goes 
further to investigate in detail the microstructure and the dynamics of a chain in the chain 
redissolution region, which gives us a thorough vision of the behavior of polyelectrolytes 
in salt solutions. Since chain length is also a relevant factor to affect the morphology of a 
chain, it will be very interesting to study the effect of finite chain length. That will be the 
topic of our future work. 
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Figure captions 

1 /2 

FIG 1 (i?^) as a function of Cg for the chains of different stiffness ka- 

FIG 2 a = {(^Rgy^^ I (i?^ — 1 as a function of Cg for the chains of different stiffness ka- 

FIG 3 q; as a function of ka at different salt concentration Cg. 

FIG 4 (-Rg)^''^ as a function of Cg for the chains of different ka- The dashed curves denote the 
prediction of the WLC model combined with the OSF and MO theories for different 
values of ka-, indicated on the left-hand side of the corresponding curve. 

FIG 5 Snapshots of simulation for a flexible chain {ka — QksT /ra(P) at (a) Cg — 1.33 x 
10-V-^ (b) C, = 8 X 10-V-^ and (c) Cg = 5.12 x IQ-^a'^. The gray bead- 
spring chain represents the polyclcctrolytc, the black spheres represent the tctravalent 
counterions, and the white spheres represent the monovalent counterions. Colons are 
not shown for the reason of clarity. 

FIG 6 Snapshots of simulation for a semiflexible chain {ka = ^ksT /radP) at (a) Cg = 1.33 x 
10-V-^ (b) Cg = 6.4 X lO-V-3, and (c) Cg = 5.12 x lO'V-^. The gray bead- 
spring chain represents the polyelectrolyte, the black spheres represent the tetravalent 
counterions, and the white spheres represent the monovalent counterions. Colons are 
not shown for the reason of clarity. 

FIG 7 Snapshots of simulation for a rigid chain {ka = 100kBT/ra(f) at Cg = 8x 10~^c7~^. The 
gray bead-spring chain represents the polyelectrolyte, the black spheres represent the 
tetravalent counterions, and the white spheres represent the monovalent counterions. 
Colons are not shown for the reason of clarity. 

FIG 8 (a) p{A), and (b) p{Rg) of the flexible chain with zero stiffness for the cases of different 
Cg- Each case is numbered in Roman number along 'Case'-axis. The value in the 
parenthesis following the Roman number denotes the salt concentration in the unit 
lO-V-3. 

FIG 9 (a) p{A), and (b) p{Rg) of the semiflexible chain with ka — SksT/racP for the cases of 
different Cg. Each case is numbered in Roman number along 'Case'-axis. The value 
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in the parenthesis following the Roman number denotes the salt concentration in the 
unit 



FIG 10 (a) p{A), and (b) p{Rg) of a rigid chain with ka — lOOkBT/racP for the cases of 
different Cg. Each case is numbered in Roman number along 'Case'-axis. The value 
in the parenthesis following the Roman number denotes the salt concentration in the 
unit 

FIG 11 Time evolution of (a) A, and (b) Rg and Re, at = 8 x lO^^cr"^ where the chain 
exhibits a hairpin structure. 

FIG 12 Time evolution of (a) A, and (b) Rg and Re, ed, Cg — 8 x 10~^(7~^ where the chain 
shows a toroidal structure. 

FIG 13 Time evolution of (a) A, and (b) Rg and Re, at Cg — 6.4 x 10~'^cr~^ where the chain 
alters between a hairpin and a toroid structure. 

FIG 14 Time evolution of (a) A, and (b) Rg and Re, at Cg = 1.28 x 10~^(T~^ where the chain 
morphology alters between a hairpin, a toroid, and an unfolded-chain structure. 

FIG 15 Schematic state diagram for our single polyelectrolytes. A double-headed arrow ap- 
peared between two states denotes that the structural transition can take place in a 
typical simulation run at that condition. 
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